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PREFACE 


The  work  on  this  program  was  divided  into  two  tasks  which 
can  be  summarized  as: 

Task  1.  A review  of  penetrator  computations  done 
by  others , and 

Task  2.  A suite  of  computations  utilizing  knowledge 
gained  in  Task  1 and  PacTech  techniques. 

The  work  on  Task  1 was  previous  reported  in,  PT-U75-0013  dated 
April  30,  1975.  This  volume  includes  the  computations,  analysis 
and  conclusions  of  Task  2. 
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1.  INTRODUCTION  AND  CONCLUSIONS 


A number  of  system  concepts  have  recently  been  proposed  which 
utilize  a wide  variety  of  delivery  and  guidance  systems  to  deliver 
a tactical  nuclear  device  to  a target  area  and  penetrate  the  ground 
surface  prior  to  detonation.  The  phenomenology  discussed  in  this 
report  is  a required  step  in  the  evaluation  of  the  loads  that  the 
vehicle  and  weapon  package  must  survive  prior  to  detonation. 

If  even  shallow  penetration  is  achieved,  prompt  radiation 
is  suppressed  as  both  a kill  and  collateral  damage  effect.  The 
alteration  of  the  remaining  primary  kill  and  collateral  damage 
effects  (air  blast  and  ground  motion)  can  be  pretested  to  a limited 
extent  with  conventional  explosives.  The  main  disadvantage  of 
penetrating  weapons  (compared  to  an  air  burst)  is  the  enhancement 
of  fallout.  A slight  reduction  in  air  blast  kill  range  is  conceded, 
but  ground  motions  are  substantially  enhanced. 

Systems  utilizing  deeper  penetrations  (DOB)  are  proposed  for 
use  when  the  target  is  sufficiently  hard  (such  as  a runway  or  buried 
and  surface  fortifications)  that  the  kill  mechanisms  of  major  inter- 
est are  the  ground  motions  (including  cratering)  and  debris  near  or 
in  the  crater.  The  fallout  hazard  can  be  substantially  reduced  by 
this  additional  burial. 

As  the  DOB  is  further  increased,  one  of  two  phenomena  occur. 

In  a medium  such  as  granite  or  basalt  which  can  support  a large 
amount  of  bulking,  a mound  will  form.  In  a weaker  material,  the 
initial  cavity  collapses;  and  a crater  with  a radius  and  volume 
approximately  equal  to  the  cavity  radius  and  volume  is  formed.  In 
either  case  the  cavity  gases  are  essentially  contained  so  there  is 
no  venting.  The  target  kill  mechanism  is  the  resulting  ground 
motion.  This  is  by  far  the  cleanest  of  the  proposed  EPW  systems. 
However,  design  feasibility  and  the  number  of  potential  targets 
have  as  yet  not  been  clearly  defined. 
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To  investigate  the  phenomenology  of  earth  penetrating  weapons 
(EPW),  the  Defense  Nuclear  Agency  (DNA),  conducted  a series  of 
experiments  at  the  Defense  Research  Establishment  Suffield  (DRES) 
Watching  Hill  Site  in  Alberta,  Canada  in  mid-July,  1974.^  Two 
vertical  firings  of  400-pound  penetrators,  five  feet  in  length, 

6.5  inches  in  diameter  were  made  at  velocities  of  300  and  500  ft/s. 
The  projectiles  were  instrumented  to  orovide  deceleration  history, 
and  stress  gauges  were  emplaced  in  the  soil  near  the  projectile's 
trajectory  to  measure  the  stress  developed  in  the  target. 

Prior  to  the  experiments,  the  Waterways  Experiment  Station 
had  collected  a great  deal  of  data  from  nearby  borings  and  excava- 
tions, and  constructed  a recommended  set  of  material  property  data. 
In  response  to  the  calculators'  request,  the  soil  was  divided  into 
the  four  most  representative  layers.  The  top  three  layers  turned 
out  to  be  eight  feet  thick.  Representative  oroperty  data  was 
supplied  by  WES  for  each  layer. This  data  was  supplied  to 
several  calculators,'  who  performed  detailed  predictive  calcu- 

lations of  the  projectile  deceleration,  penetration  depth,  and 
stresses  on  the  projectile  and  in  the  soil  for  the  case  of  the  test 
with  the  nominal  500  ft/s  impact  velocity.  The  calculations  using 
empirical  and  analytical  techniques  are  reviewed  and  compared  in 
Reference  13.  Those  using  finite  difference  techniques  have  been 
reviewed  and  compared  in  Reference  14. 

From  a post-shot  excavation  made  at  the  shot  site,  WES 
discovered  that  the  top  3 to  5 feet  of  the  firing  site  had  been 
regraded  and  was  composed  of  fill  material,  rather  than  of  the  top 

(13) 

layer  (Layer  1)  material  which  they  had  supposed.  ' In  addition, 
the  deceleration  measurements  disclosed  that  the  material  near  a 
depth  of  12  feet  had  characteristics  which  created  especially  high 
decelerations.  ^ 
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The  computations  of  this  study  consider  the  materials  of 
Layer  1,  and  of  the  fill  layer  discovered  postshot.  The  sequence 
of  computations  was  laid  out  so  as  to  permit  comparisons  between 
cases  in  which  only  one  parameter  had  been  altered,  which  provides 
an  estimate  of  the  relative  importance  of  each  parameter.  The  suite 
of  computations  is  discussed  in  Section  2. 

The  inputs  to  the  computations  are  discussed  in  Section  3, 
including  the  projectile  geometry  and  mass,  and  the  material  pro- 
perties used  for  the  fill  and  Layer  1.  The  computational  results 
are  discussed  in  Section  4. 

T'  jits  indicate  that  the  fill  material  produces  less 
dece1  *:han  the  material  of  Layer  1,  so  that  the  sequence  of 

lay  the  projectile  passes,  viz.  fill  to  Layer  1,  provides 

an  . i ng  sequence  of  decelerations.  This  provides  a partial 

explanation  for  the  apparent  increase  of  deceleration  with  depth, 
which  was  initially  a puzzling  aspect  of  the  experimental  measure- 
ments. To  complete  the  explanation  of  this  phenomenon  the  proper- 
ties of  Layer  2 should  be  more  carefully  investigated.  If,  as 
expected,  this  results  in  higher  computed  decelerations  than  in 
Layer  1,  a complete  interpretation  of  this  phenomenon  would  result. 

All  of  our  computations  indicated  that  a "steady-state"  flow 
field  had  arisen  about  the  projectile  after  about  18  inches  of 
penetration,  after  which  the  deceleration  depends  only  on  the  pro- 
jectile velocity,  which  changes  relatively  slowly  with  depth.  Thus, 
it  is  necessary  to  compute  only  about  24  inches  of  penetration  in 
a given  material  to  establish  the  equilibrium  value  of  deceleration 
at  a given  velocity.  Conducting  a sequence  of  such  computations  at 
various  velocities  would  therefore  provide  a means  to  compute  the 
projectile  deceleration  and  loadings  as  a function  of  velocity  for  a 
given  medium. 

The  zoning  study  indicated  that  three  centimeter  zones  are  not 
small  enough  to  permit  an  accurate  computation  of  the  deceleration 
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and  loading  for  projectiles  with  the  blunted  nose  tips  considered 
in  this  study.  The  inadequacy  becomes  worse  at  higher  impact 
velocities.  The  zoning  study  should  be  continued  to  much  smaller 
zone  sizes  to  determine  the  required  zoning  fineness. 

In  Cases  2 through  5,  the  calculations  utilized  fits  to  the 
data  furnished  by  WES,  and  did  not  include  frictional  forces  between 
the  soil  and  projectile.  Using  the  material  property  data  in  this 
way  will  ultimately  lead  to  a definite  computed  deceleration  when 
sufficiently  fine  zones  are  used.  Although  these  final  values  of 
deceleration  are  unknown  at  present,  the  values  generated  in  this 
study  suggest  that  little  or  no  friction  will  have  to  be  added  to 
bring  the  computed  decelerations  into  correspondence  with  the  measured 
values.  This  conclusionof  course  applies  only  to  the  materials 
treated  in  this  study,  namely  fill  and  Layer  1 material.  In  addition, 
it  assumes  that  the  measurements  of  deceleration  were  accurately 
conducted,  and  it  should  be  pointed  out  that  the  measurements  contain 
inexplicable  oscillations  which  are  still  under  investigation. 

The  calculations  described  herein,  a first  exploratory  cut, 
were  funded  at  the  level  of  a sixth  of  a man-year.  They  reveal  that 
finer  zoning  should  be  used  to  obtain  quantitative  improvements 
in  the  calculated  stress  on,  and  deceleration  of,  a blunted-tipped 
projectile.  The  computational  techniques  employed  in  this  study  have 
proved  to  be  so  efficient  that  it  is  easily  feasible  to  use  much 
finer  zoning.  All  cases  studied  used  1500  cells  and  required 
expenditures  of  between  306  and  625  seconds  of  central  processing 
(CPU)  time  on  the  CDC  7600  to  compute  four  feet  of  penetration 
(5.1  to  10.4  minutes).  Cells  near  the  penetrator  were  integrated 
for  4000  to  8000  cycles,  whereas  those  far  from  the  penetrator 
required  only  300  to  600  cycles.  The  code  contains  an  automatic 
subcycling  routine  which  selects  the  optimum  region  to  be  sub- 
cycled to  advance  the  solution  farthest  in  real  time  for  a given 
expenditure  of  computer  time.  In  addition,  the  use  of  automatic 
rezoning  permitted  the  problems  to  be  run  to  completion  without 
manually  rezoning  the  comDuting  grid. 
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2.  THE  CALCULATIONS 


Five  calculations  were  conducted  in  this  study,  all  using 
the  vehicle  geometry  described  in  Section  3.1. 

Case  1 . 


The  soil  used  in  this  calculation  had  the  properties 
originally  specified  by  WES  for  Layer  1 of  the  Watching  Hill  Site. 
The  impact  velocity  was  500  ft/sec.  Friction  was  included  in  the 
form 


Jo  ) 

| t ! = Min  {0.6o  , 0.6  4^1 
n v j 


(2.1) 


t being  applied  so  as  to  oppose  the  differential  velocity  between 
the  soil  and  projectile.  Near  the  penetrator,  the  computing  zones 
were  initially  square  in  cross-section  having  sides  of  4.1275  cm. 
This  computation  modeled  as  closely  as  possible  that  conducted  by 
CRT. (10)  Toward  this  end  the  CRT  material  property  fit  was  used. 

Case  2. 

The  specifications  for  this  calculation  are  identical  with 
those  of  Case  1 with  a single  exception.  No  friction  was  included 
in  the  model. 

Case  3. 

This  calculation  is  identical  with  Case  2 in  every  respect 

except  that  the  soil  had  the  characteristics  of  the  fill  material 

which  was  discovered  at  the  top  of  the  Watching  Hill  Site  after  the 
f 1 31 

shot.  The  material  properties  of  this  material,  and  the  fit 

used  in  the  calculation  are  discussed  in  Section  3.2. 

Case  4. 

This  calculation  is  identical  in  all  respects  with  Case  3, 
except  that  an  impact  velocity  of  3500  ft/sec  was  used. 
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TABLE  2.1 


SUMMARY  OF  COMPUTATIONS 

Case 

1 

2 

3 

4 

5 

Friction 

yes 

no 

no 

no 

no 

Soil 

Layer  1 

Layer  1 

Fill 

Fill 

Layer  1 

Velocity  (ft/s) 

500 

500 

500 

3500 

500 

Zone  size  (cm) 

4.1275 

4.1275 

4.1275 

4.1275 

3 

Decel  (g's) 

191  .2 

89.5 

59.5 

1982 

81.9 

Approximate  measured 

deceleration  (g's)* 

75 

75 

20-90** 

NA 

75 

*From  Figure  5.1-1  (See  Reference  13) 

**There  is  an  unexplained  oscillation  in  the  data  from  the  500  ft/s 
test  v/h i 1 e the  projectile  is  traversing  the  fill.  The  estimated 
average  value  is  50  g's. 


3.  INPUTS  TO  THE  COMPUTATIONS 


3.1  VEHICLE  GEOMETRY 

Figure  3.1-1  shows  a drawing  of  the  actual  vehicle  together 
with  the  manufacturing  tolerances.^  The  cylindrical  afterbody 
has  a diameter  of  6.5  in.,  and  a length  of  41.524  ± .03  in. 

The  nose  is  ogival  in  character,  being  specifically  a 9.25  caliber 
radius  head  (CRH)  tangent.  This  means  the  outline  of  the  nose  is 
generated  by  a circle  of  radius  9.25  projectile  diameters  or 
60.125  in.  The  circle  is  tangent  to  the  afterbody  at  the  point 
where  it  intersects  it.  As  the  figure  shows,  this  implies  that 
the  circle  center  is  at  the  same  station  as  the  forward  part  of  the 
afterbody,  and  is  offset  from  the  vehicle  centerline  by  8.75  calibers 
or  56.875  in.  If  the  nose  were  a complete  ogive,  it  would  have  a 
length  of  19.5  in.  However,  the  nose  is  blunted  by  the  addition 
of  the  conical  tip  0.5  in.  high  having  a 45°  half  angle.  Milling 
off  the  conical  tip  reduces  the  length  of  the  ogive  by  2.602  cm, 
or  1.024  in.,  leading  to  an  actual  nose  length  of  18.476  in.  The 
length  of  the  afterbody  plus  the  length  of  the  blunted  nose  thus 
totals  a nominal  60  in.  As  seen  from  the  figure,  the  manufacturing 
specifications  are  that  vehicle  be  60.00  + .06  in.  in  length. 

Figure  3.1-2  shows  two  detailed  drawings  of  the  nose  geometry. 
Figure  3 . 1 -2a  is  a full-scale  drawing  of  the  first  20  cm  of  the 
nose.  Figure  3 . 1 - 2b  is  full  scale  in  the  radial  direction,  but 
the  Y-scale  (distance  from  tip)  is  compressed  by  a factor  of  2. 

In  the  calculations  to  be  presented,  the  projectile  is  treated  as 
a rigid  body  so  that  it  always  preserves  the  dimensions  and 
configuration  shown  in  the  figure. 

The  weight  of  the  vehicle  is  a nominal  400  lbs.,  and  the 
equivalent  value  of  181,436  grams  is  used  in  the  calculation. 

In  the  calculation,  the  afterbody  was  terminated  flush  at  a distance 
of  60  in.  from  the  nose  tip.  This  neglects  the  screw-in  insert 


10 


Figure  3.1-2b.  Nose  tip  geometry  (Y:  half  scale,  R:  full  scale) 


shown  in  Figure  3.1-1.  However,  it  was  observed  that  the  soil 
separated  from  the  vehicle  in  the  computation  and  contacted  no 
portion  of  the  afterbody,  so  that  the  flush  termination  did  not 
influence  the  computed  deceleration. 
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3.2  MATERIAL  PROPERTIES  OF  THE  FILL  LAYER 


3.2.1  The  Data 

(151 

Based  on  experimental  data  provided  by  WES,  ' an  analytical 
fit  to  the  material  properties  of  the  fill  layer  was  conducted 
for  use  in  these  computations.  The  data  transmitted  by  WES  was  a 
series  of  recommended  curves,  based  on  the  results  of  experiments 
made  on  several  samples.  The  final  recommendations  were  constructed 
using  averages  of  the  test  results,  combined  with  judgment  as  to 
how  the  test  data  related  to  the  material  properties  in  situ. 

Figures  3.2-1  through  3.2-4  contain  the  WES  recommendations, 
concerning  the  dynamic  constitutive  properties  of  the  backfill 
material  encountered  near  the  ground  surface  in  DNA  Penetration 
Test  06,  which  was  conducted  near  the  PRAIRIE  FLAT  HEST  TEST 
Site  at  the  Watching  Hill  Blast  Range  of  the  Defense  Research 
Establ ishment,  Suffield  (ORES),  Alberta,  Canada. 

The  data  transmitted  by  WES  consisted  of  the  following: 

1.  A stress-strain  loadinq  curve  in  uniaxial  strain 

8 2 

up  to  a peak  stress  of  a = 3.41  x 10  dyne/cm 
(4946  psi ) . 

2.  Stress-strain  curves  for  unloading  in  uniaxial 
strain,  beginning  at  various  points  on  the  UX  loading 
curve.  Three  such  curves  were  provided  which  begin 

on  the  loading  curve  at  a = 400,  1000,  and  4946  psi. 

3.  Stress  difference  (a  - a ) vs.  mean  stress  curves 

z r 

for  loading  in  uniaxial  strain  up  to  a mean  stress 
of  425  psi . 

4.  Curves  of  stress  difference  vs.  mean  stress  for 
unloading  in  uniaxial  strain.  Two  such  curves  were 
provided  corresponding  to  unloading  from  the  uniaxial 
strain  loading  curve  at  o = 400  and  1000  psi. 
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5. 


A failure  envelope  couched  in  terms  of  principal 
stress  difference  as  a function  of  mean  stress  as 
deduced  from  triaxial  strain  tests.  Data  was 
presented  up  ; a mean  stress  of  1650  psi. 

The  curves  in  Fic.u“  ■ 3.2-1  through  3.2-4  are  in  terms  of 
live  stresses  (i.e.,  the  u situ  state  of  stress  corresponds  to 
the  origin  of  coordinates  in  the  plot).  The  strains  in  these 
figures  are  expressed  as  engineering  strains.  Both  stress  and 
strain  are  considered  positive  in  compression. 

The  depth  of  the  backfill  in  the  immediate  vicinity  of 

Test  06  varies  from  3 to  5 feet.  A backfill  layer  thickness  of 

(15) 

4 feet  was  recommended  by  WES  for  use  in  these  calculations 
simulating  the  penetration  event.  Below  the  backfill,  the 
recommended  pretest  Layer  1 properties  still  appear  to  be  appropriate 
for  the  4 to  8 foot  depth  range. 

Laboratory  test  results  from  undisturbed  samples  obtained 
from  the  backfill  zone  in  the  post-test  field  investigation  were 
used  as  the  basis  for  developing  the  properties  contained  in 
Figures  3.2-1  through  3.2-4.  These  tests  consisted  of  uncon- 
sol idated-undrained  uniaxial  strain  (UX)  and  triaxial  compression 
(TX)  tests.  The  test  specimens  exhibited  a rather  broad  range 
of  scatter  in  their  water  contents  and  densities.  This  is  not 
unusual  for  a material  of  this  type  placed  without  any  special 
controls.  It  is  WES's  judgment  that  the  backfill  should  be 
treated  as  a uniform  layer  with  the  following  representative 
composition  properties: 

Water  Content  (weight  of  H^O  per  total  wet  weight)  = 18" 

Wet  Unit  Weight  = 97  lb/ft^  = 1.55449  g/cm^ 

The  triaxial  tests  conducted  by  WES  contain  uncertainties 
due  to  end-cap  friction.  The  axial  and  radial  strains,  e , and  of 
are  respectively  determined  from  the  overall  shortening  of  the 
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specimen,  and  the  radial  deformation  measurement  at  the  specimen  mid- 
point. To  accommodate  the  end-cap  uncertainty,  WES  computed  the 
volumetric  comoaction,  Av/vQ  in  two  different  ways.  The  first  way 
assumes  that  end-cap  friction  is  negligible,  and  that  the  cylin- 
drical specimen  remains  cylindrical.  This  leads  to  the  equation 


Av  _ ,0 

— = f.  + 2r 

v a r 
o 


Zr  r + e2  (c  - 1)  + 
r a r a 


(3.2-1) 


The  second  way  assumes  that  high  end-cap  friction  exists  and  that 
the  deformed  specimen  takes  the  shape  of  a double  conic,  pinched 
at  the  midpoint.  The  volumetric  compaction  in  this  case  is 
calculated  from 


t.  r 

r a 


(r  - 1) 


(3.2-2) 


When  applied  to  a given  set  of  data  for  the  shear  phase  of  a triaxial 
test,  the  first  equation  will  tend  to  minimize  apparent  compaction 
and  maximize  apparent  dilation.  The  latter  formula  has  the  opposite 
effect.  Examining  the  triaxial  data  with  these  factors  in  mind, 

WES  concluded  that  the  backfill  materials  tested  definitely  do  not 
exhibit  any  significant  dilation  during  straining  in  shear.  They 
either  compact  or  remain  at  approximately  constant  volume  depending 
on  which  deformed  shape  is  assumed.  In  view  of  this,  an  elastic, 
ideally-plastic  (EIP)  type  of  constitutive  model  with  a non- 
associated  flow  rule  is  believed  to  be  physically  more  appropriate 
than  one  with  an  associated  flow  rule. 

In  constructing  the  recommended  curves  shown  as  Figures 
3.2-1  through  3.2-4,  WES  intentionally  made  the  curves  stiffer 
than  any  of  the  data  actually  taken.  This  was  done  to  account 
for  differences  in  initial  density,  loading  rate,  and  possible 
sample  disturbance.  These  changes  manifest  some  of  the  judgments 
which  WES  had  to  apply  in  formulating  these  recommended  curves. 
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RECOMMENDED  TX  FAILURE  SURFACE 
RECOMMENDED  I IX  STRESS  PATH 


d TX  failure 


1 


isd  *(y.  - 2 i)  93uaj3^>ig  ssaj^s  ledpuug 


Figure  3.2-4.  Recommendation  of  triaxial  test  failure  surface  for  the  backfill. 


The  recommended  triaxial  strain-based  failure  surface  for 
the  backfill  material  (Figures  3.2-3  and  3.2-4)  represents  a 
considerable  reduction  in  shear  strength  from  that  of  the  pretest 
Layer  1 material.  This  curve  was  intentionally  set  above  the 
available  triaxial  test  data  points  because  the  test  specimens 
had  water  contents  higher  than  those  judged  to  be  representative 
of  the  entire  stratum. 

The  stress  paths  for  the  uniaxial  strain  test  shown  in  Figure 

3.2- 4  correspond  to  the  uniaxial  stress-strain  curves  shown  in  Figure 

3.2- 1.  No  stress  path  data  were  available  on  which  to  base  the  re- 
commended curves,  however,  they  are  qualitatively  similar  to  those 
recommended  for  the  pretest  Layer  1 material  and  the  initial  slope 
correlates  with  the  expected  values  calculated  from  the  initial 
moduli  data  extractable  from  the  individual  test  results. 
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3.2.2  The  Fit 


The  experimental  data  provided  by  WES  are  not  sufficient  to 
determine  the  equation  of  state  of  the  material  completely.  The 
penetration  process  causes  the  soil  to  traverse  different  stress- 
strain  paths  than  those  along  which  data  were  provided.  It  is 
therefore  necessary  to  make  additional  assumptions  as  to  the  con- 
stitutive relations  in  order  to  generate  a material  property  model 
which  applies  to  the  region  of  stress-strain  space  encountered 
during  penetration.  In  the  absence  of  further  data,  these 
assumptions  must  be  arbitrary  to  a great  extent.  They  are 
restricted  only  by  the  requirements  that  the  constitutive  relations 
must  satisfy  the  criteria  of  uniqueness  and  thermodynamic  con- 
sistency. However,  these  requirements  are  very  broad,  as  they 
must  be  to  encompass  every  material  in  the  world! 

Specifically,  the  additional  major  assumptions  were  made 
that  the  loading  and  unloading  bulk  moduli  are  functions  only 
of  pressure,  and  that  the  loading  and  unloading  rigidity  moduli 
are  functions  only  of  v4T^'.  Using  these  assumptions,  the  loading 
and  unloading  curves  in  uniaxial  strain  are  considered  to  be 
reference  curves  which  provide  data  to  fit  the  loading  and  unload- 
ing modul i . 

The  reference  curves  used  to  find  the  loading  moduli  are 
the  UX  o-c  loading  curves,  and  UX  Y-P  loading  curves.* 

As  shown  in  Figure  3.2-2,  the  o-c  loading  curve  provides  data 
up  to  oz  = 4946  psi  (3.41  x 10  dynes/cnr).  The  curve  was 
arbitrarily  extrapolated  to  higher  values  of  stress  by  assuming 
continuity  of  the  first  derivative,  and  constancy  of  the  second 
derivative.  The  UX  Y-P  loading  curve  is  shown  in  Figure  3.2-3. 


*To  simplify  reference  to  the  data,  the  term,  "UX  o-c  curves", 
refers  to  the  stress-strain  paths  taken  in  uniaxial  strain 
(Figures  3.2-1  and  3.2-2),  while  "UX  Y-P"  refers  to  the  stress 
difference  vs.  mean  stress  plots  (Figure  3.2-4). 
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It  extends  only  up  to  a mean  stress  of  425  bars,  so  that  it  had 
to  be  arbitrarily  augmented  above  that  mean  stress.  At  a mean 
stress  of  425  bars,  it  falls  about  5 psi  beneath  the  failure 
surface.  It  was  assumed  that  it  remained  offset  by  5 psi  from 
the  failure  surface  which  is  given  up  to  a mean  stress  of  1650  psi 
in  Figure  3.2-4.  Both  the  failure  surface  and  UX  Y-P  loading 
curve  were  further  augmented  by  assuming  that  the  failure  surface 
approaches  the  value  of  250  psi  asymptotically  at  high  mean 
stress.  The  offset  of  5 psi  between  the  UX  Y-P  loading  curve 
and  the  failure  surface  was  preserved  throughout  the  entire  region 
of  higher  mean  stress. 

The  reference  curves  used  to  determine  the  unloading  moduli 
are  1)  the  unloading  UX  o-c  curve  shown  in  Figure  3.2-2,  which  is 
tangent  to  the  loading  curve  at  4946  psi,  and  2)  the  UX  Y-P  curve 
shown  in  Figure  3.2-4  which  encounters  the  UX  loading  curve  at 
cr  = 1000  psi.  It  was  assumed  that  curve  1)  is  the  limiting 
unloading  curve,  and  that  it  is  tangent  to  the  UX  loading  curve  at 
o = 4946  psi,  (Point  T in  Figure  3.2-2).  If  UX  loading  to  a 
higher  stress  occurs,  the  material  will  unload  back  along  the 
loading  curve  until  Point  T is  reached,  and  it  will  then  follow 
the  limiting  unloading  curve.  Point  T is  the  point  at  which  all 
voids  are  irreversibly  crushed  out. 

Data  on  the  Y-P  unloading  curve  was  provided  only  up  to  a 
mean  stress  of  425  psi,  and  had  to  be  arbitrarily  augmented  at 
higher  mean  stress.  This  was  done  in  a series  of  steps.  First, 
the  augmented  UX  Y-P  loading  curve  was  used  to  determine  the 
stress  difference  which  occurs  in  UX  strain  when  a = 1000  psi.  j 

The  reference  unloading  curve  was  then  extrapolated  from  a mean 
stress  of  425  psi  to  this  point,  preserving  the  continuity  of 
its  first  derivative  at  P = 425  psi. 

Having  specified  these  curves  up  to  a mean  stress  of  1000  psi, 
it  is  possible  to  determine  the  rigidity  modulus  up  to  the  stress 
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difference  which  occurs  at  1000  psi.  There  remains  the  problem 
of  defining  it  at  higher  values  of  the  stress  difference,  or  more 
generally,  higher  values  of  /J^ ' . This  was  done  in  the  following 
way.  It  was  assumed  that  at  a mean  stress  of  4946  psi,  the  voids 
are  all  crushed  out,  and  the  loading  and  unloading  curves  should 
be  tangent  at  that  point  and  identical  at  higher  mean  stresses. 
Thus,  both  G and  its  derivative  with  respect  to  are  defined 
at  Point  T.  The  previous  procedure  defined  both  G and  its  first 
derivative  with  respect  to  /J^'  at  a mean  stress  of  1000  psi. 
Therefore,  a simple  polynomial  fit  was  used  to  span  the  gap, 
preserving  the  continuity  of  G and  its  first  derivative  in  the 
interval  up  to  the  point  at  which  the  voids  are  crushed  out. 

To  construct  our  material  model,  we  first  fit  o|"  (t:  ) 

U z z 

and  o2  (rz)  which  are  the  loading  and  unloading  curves  in  uniaxial 

strain.  In  this  notation,  it  is  assumed  that  the  strain  occurs  in 

the  Z direction,  cz  is  the  strain  in  that  direction,  and  o 

the  associated  stress.  We  next  use  the  uniaxial  strain  loading 

and  unloading  curves  to  fit  Y^(P)  and  Y^(P),  where  Y is  defined 

through 


v2  = 2 

2 sij  si j 

so  that  in  the  case  of  uniaxial  strain. 


Y = o-o 

1 z r 1 


(3.2-3) 


(3.2-4) 


On  loading,  a is  always  greater  than  or,  and  on  unloading  we 
conduct  the  fit  only  to  the  point  where  Y first  becomes  zero. 
Thus,  a is  greater  than  o throughout  this  entire  loading- 
unloading cycle.  For  this  loading  and  unloading  region,  we  then 
have  the  identity 


z 


(3.2-5) 


This  identity,  and  the  fits  of  °z(!z)  an^  Y(P),  permit  us  to 


solve  for  P as  a function  of  c 


Since  t:  depends  only  on 
z J 


the  volumetric  strain  in  these  uniaxial  strain  tests,  this  pro- 
vides us  with  our  pressure-volume  relation. 


Because  the  engineering  strain,  as  defined  by 


v - v 
o 


(3.2-6) 


is  used,  the  bulk  modulus  is  prescribed  by 


n do  1 - y 

K = - v -P  = — 2 L_- 

v :»v  d>  _ , . 2 dY 


(3.2-7) 


2 1 + 5 dP 


To  obtain  the  rigidity  modulus  from  these  fits,  we  start 
with  the  definition 


i / v , 4G , d v 

do  = - (K  + v -) 
z 3 v 


dP  = - K 


dv 


(3.2-8) 


dv 


where  — - is  the  true  incremental  strain  as  used  in  our  computations. 


This  leads  to  the  definition 


1 AD  d'1  1 (1  " r d'1 

r _ 3 ,dP  z\  _ 1 ' z dP  _ z 
b " 4 v dv  ‘ d v"  ' 2 ; . 2 dY  dt  ' 

' 3 d P z 


(3.2-9) 


To  obtain  the  loading  rigidity  modulus,  the  loading  fits  are  used, 
and  the  unloading  modulus  utilizes  the  unloading  fits. 


This  approach  permits  us  to  conduct  fits  of  the  data  itself, 
and  then  deduce  the  moduli  from  the  fit.  In  conducting  the  fit 
it  is  necessary  to  insure  the  continuity  of  both  the  function 
and  its  first  derivative.  This  is  accomplished  by  breaking  the 
independent  variable  into  regions,  conducting  fits  in  each  region 


wherein  the  function  and  its  first  derivative  are  piece-wise 
continuous,  and  matching  the  function  and  its  first  derivatives 
at  the  intersection.  Quadratic  functions  are  especially  suitable 
for  this  purpose,  because  they  can  easily  be  inverted. 

The  approach  has  the  virtue  that  if  a loading-unloading  calculation 
in  UX  strain  is  conducted,  the  experimental  data  will  be  precisely 
reproduced.  Whereas  the  resulting  moduli  are  continuous,  their 
derivatives  may  be  discontinuous.  However,  the  discontinuties 
in  the  moduli  derivatives  leads  to  no  difficulties. 

The  data  provided  by  WES  for  the  failure  locus  is 
(o  - o ) as  a function  of  the  mean  stress.  To  accommodate  this 
description,  we  define  a variable,  Y through  the  relation 


s..  s..  = | Y2 
i.l  ij  3 


(3.2-11) 


where  s.  . are  the  physical  components  of  the  stress  deviator. 
Because  of  the  definition 


• V * •T's-.s.T 

£ 2 ij  i.) 


Y and  are  related  through 


Y = /3  J2‘ 


(3.2-12) 


(3.2-13) 


For  loadings  in  which  the  two  minor  principal  stresses  are  equal, 
Y becomes  just  the  difference  between  the  greatest  and  the  least 
principal  stresses.  Thus,  for  the  case  of  uniaxial  strain,  or 
triaxial  loading,  we  have  the  relation 


Y = o -a 
i z r 1 


(3.2-14) 


We  fit  Y as  a function  of  mean  stress  by  means  of  a series 
of  quadratic  expressions,  maintaining  continuity  of  the  function 


and  its  first  derivative.  The  fit  has  the  general  form 


YF  = BQ  + B.jP  + B2P2  (3.2-15) 

The  coefficients  used  in  the  fit  are  given  in  Table  3.2-1.  The 
I-th  set  of  coefficients  is  used  when 

PT ( I - 1)  < P < PT( I ) , (3.2-16) 

The  set  corresponding  to  I = 1 is  used  when  P < PT  (1). 

In  the  case  of  unconfined  compression,  the  two  minor  stresses  are 
zero,  and  Y and  P are  related  through 

oz  - or  = 3P  (3.2-17) 

If  o is  greater  than  a , this  has  the  form 


YF  = 3P  (3.2-13) 

The  unconfined  yield  stress  is  prescribed  by  the  intersection  of 
the  line  defined  by  Equation  (3.2-18)  and  the  parabolic  section 
defined  by  Equation  (3.2-5)  with  I = 2.  The  unconfined  yield 
stress  and  corresponding  mean  stress  found  in  this  way  are 


Y = 1.5712  x 106  dynes/cm2  = 22.79  psi 
Py  = 5.376  x 10®  dynes/cm2  = 7.60  psi 


(3.2-19) 


For  a mean  stress  more  tensile  than  10  osi  = 6.894757  x 10 
2 

dynes/cm  Y is  set  uniformly  to  zero. 
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TABLE  3.2-1 


CONSTANTS  FOR  THE  FIT  T 

0 THE  FAILURE  LOCUS 

I 

P(I) 

, fyn  , 

(dyne/cm^) 

Bl(I) 

dimensionless 

MI) 

(cnr/dyne) 

1 

-6.89475727+05 

0.00000000 

0.00000000 

0.00000000 

2 

5.17106800+06 

9.12224844+05 

1.28626876+00 

-5.33761018-08 

3 

1 . 20658252-07 

1.35345691+06 

1.11561462+00 

-3.68752424-08 

4 

2.06842720+07 

5.58223081+06 

4.14663989-01 

-7.82830145-09 

5 

2.75790292+07 

1.11860145+07 

-1.27176087-01 

5.26957500-09 

6 

5.51580585+07 

5.11682356+06 

3.12954776-01 

-2.70987310-09 

7 

4.48159224+08 

1.33072241+07 

1 .59754809-02 

-1  .77977936-11 

8 

9.99999998+19 

1.68818520+07 

2.29902400-05 

0.00000000 

5.17+6  denotes  5.17x1 0^ , etc. 


28 


3 . 3 MATERIAL  PROPERTIES  OF  LAYER  7 

Since  Case  1 was  supposed  to  reproduce  the  calculation 
previously  conducted  by  CRT  for  the  Layer  1 as  closely  as  possible, 
the  CRT  fit  to  the  material  properties^)  was  used.  To  minimize 
the  variations  between  Case  1 and  Cases  2 and  5,  the  same  fit 
was  used  in  those  calculations  also. 

As  discussed  more  fully  in  Reference  14,  the  data  on 
which  this  fit  was  based  was  supplied  by  WES.  Figures  3.3-1 
through  3.3-6  summarize  the  WES  data,  and  compare  it  with  the 
CRT  fit. 


Vertical  Stress 


Figure  3.3-2.  Uniaxial  Strain  Load-Unload  Paths  for  Layer  1 
(az  to  14  KSI ) 


3 


Failure  Surface 


Figure  3.3-3.  Uniaxial  Strain  Stress  Paths  and  Failure  Surface  for  Layer  1 (P  to  400  psi) 


re  Surface  in  High-Pressure  Region  for  Layer 


pressu 


Volumetric  Strain  u 


Figure  3.3-5.  Hydrostatic  Load-Unload  Paths  for  Layer  1 
(P  to  1000  psi) 


4.  CALCULATION  RESULTS 

4. 1 Deceleration 

The  computated  decelerations  for  the  cases  having  impact 
velocities  of  500  ft/s  are  summarized  in  Figures  4.1-1  and 

4.1- 2.  It  was  found  that  the  velocity-time  plots  of  Figure 

4.1- 1  were  the  most  accurate  means  of  determining  the  deceleration 
as  a function  of  time.  The  computed  deceleration  curves  are 
qualitatively  similar  for  all  four  cases  shown.  After  initial 
impact,  the  deceleration  increases  smoothly  as  more  and  more  of 
the  nose  is  buried.  At  the  impact  velocity  of  500  ft/s,  the 
projectile  penetrates  about  six  inches  per  ms,  so  that  it 
requires  about  three  milliseconds  to  bury  the  nose.  It  may  be 
seen  from  the  Figure  that  a limiting  value  of  deceleration  is 
reached  by  about  three  ms,  after  which  the  velocity-time 

plot  is  essentially  a straight  line. 

It  might  seem  surprising  that  the  limiting  deceleration 
is  attained  so  rapidly.  The  reason  is  that  the  forces  which 
produce  the  deceleration  are  concentrated  near  the  nose  tip.  By 
the  time  the  ogival  portion  has  been  buried,  a "steady-state" 
flow  field  has  been  established  in  the  soil  near  the  nose  tip. 

At  that  point  in  time,  the  forces  near  the  tip  become  "constant". 
The  quotation  marks  have  been  used  because  the  flow  is  "steady- 
state"  and  the  forces  are  "constant"  on  the  time  scale  of  a few 
milliseconds  which  we  are  examining  here.  The  flow  field  and 
velocity  field  are  velocity  dependent,  and  will  change  as  the 
projectile's  velocity  decreases.  In  the  cases  we  consider  here, 
the  velocity  changes  on  the  order  of  a percent  per  millisecond, 
and  the  deceleration  and  forces  experience  changes  at  a percentage 
rate  of  the  same  order  after  "steady  state"  is  attained. 

This  assertion  is  supported  by  the  linearity  of  the  velocity-time 
plots  after  three  milliseconds. 
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VELOCITY  (10  ‘M/S) 
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Figure  4.1-1. 


TIME  (10" 3 S) 

Velocity-time  plots  for  EPW  penetrations. 
Impact  velocity  is  500  ft/s  = 152.4  m/s. 
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VELOCITY  (10  " M/S) 
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Case  2 


Layer  1 , no 
friction  ■ 
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Layer  1 , with 


friction 


j Case  5 is  same  as  Case  !2, 
! with  finer  zoning  usedj 


Case  5 
Small  zones 


V. 
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DEPTH  (M) 

Figure  4.1-2.  Velocity-depth  cross  plots  for  EPW  penetrations. 
Impact  velocity  is  500  ft/s  = 152,4  m/s. 


Deceleration  ( g * s ) 


Case 

Curve  through 
Velocity 
Peaks 

Curve  through 
Vel oci ty 
Valleys 

Average 

1 

191 .04 

191.44 

191 .24 

2 

89.79 

89.81 

89.80 

3 

59.46 

59.52 

59.49 

*4 

1982.00 

5 

81.95 

1 

81  .84 

81  .90 

*The  curves  for  Case  4 are  more  irregular,  and  only  the  average 
value  is  given.  (See  Figure  4.3-3a.) 

The  peaks  and  valleys  arise  because  of  the  finite  size 
of  the  cells  used  in  the  computing  grid.  When  the  tip  of  the 
projectile  first  encounters  a cell,  the  cell  corner  is  given 
a sudden  velocity  increment  in  one  computing  cycle.  The  velocity 
increment  is  "felt"  by  the  entire  cell  mass.  The  sudden  change 
in  cell  momentum  arising  thereby  produces  a corresponding  change 
in  projectile  momentum,  and  produces  a valley.  It  would  be  more 


accurate  to  impart  the  sudden  velocity  increment  only  to  the  mass 
of  the  cell  spanned  by  the  wave  arising  from  the  collision,  as 
would  occur  in  the  case  of  infinitely  fine  zoning.  It  is 
recognized  that  this  phenomenon  tends  to  overestimate  the  computed 
deceleration. 

The  amount  of  the  overestimate  can  be  determined  by  a 
zoning  study  in  which  the  zone  sizes  are  successively  reduced  in 
a series  of  calculations.  One  would  expect  to  find  a limiting 
zone  size  beyond  which  the  computed  deceleration  did  not  vary 
significantly.  A step  in  this  direction  was  taken  in  Cases  2 
and  5,  which  differ  only  with  respect  to  the  zone  size  used. 

Cases  2 and  5 utilized  initially  square  zones  having  sides  of 
4.1275  and  3 cm,  respectively.  The  masses  of  the  cells  initially 
on  axis  were  343.5  and  131.8  gm.  The  computed  decelerations 
were  89.8  and  81.9  g's. 

Although  the  decrease  in  computed  decleration  appears  to 
be  modest,  the  finer  zoning  was  found  to  produce  a striking  change 
in  the  distribution  of  tractions*  on  the  projectile,  as  will  be 
discussed  more  fully  in  Section  4.2.  In  view  of  this,  one  cannot 
be  certain  that  the  10%  change  in  computed  deceleration  indicates 
that  the  zoning  is  close  to  that  required  to  compute  the  correct 
deceleration.  It  is  important  to  conduct  computations  with  much 
finer  zoning  to  determine  both  the  correctly-computed  deceleration 
and  tractions*  on  the  vehicle. 


★ 

Surface  traction  is  a vector  whose  magnitude  is  the  force  per  unit 
area.  The  traction  provides  the  proper  boundary  condition  for  the 
stress  field  in  the  vehicle,  which  is  T-*  = a1Jnj,  where  T1  is  the 
traction,  o^J  is  the  stress  in  the  vehicle  evaluated  at  the  surface, 
and  nj  is  the  unit  outward  vector  normal  to  the  surface. 
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4 . 2 Stress  on  the  Projectile 

As  discussed  in  the  previous  section,  the  coarse  computing 
cell  size  leads  to  irregularities  in  the  deceleration.  The 
same  phenomenon  leads  to  irregularities  in  the  stress  and  velocity 
fields  in  the  surrounding  soil.  As  a result,  it  can  be  somewhat 
misleading  to  examine  a snapshot  of  the  stress  and  velocity  fields 
at  a particular  time.  It  is  much  more  accurate  to  conduct 
suitable  time  averages  to  determine  the  true  values  of  these 
variables.  Such  time  averages  were  conducted  to  find  the  stress 
on  the  projectile. 

The  averaging  process  is  made  more  complex  by  the  fact 
that  the  Lagrangian  computing  cells  move  along  the  Drojectile's 
surface.  To  accommodate  this,  the  surface  was  divided  into 
segments  one-half  inch  (1.27  cm)  in  height,  and  the  surface 
forces  on  each  segment  were  monitored  in  time.  The  stations 
of  the  segment  boundaries  are  tabulated  in  Table  4.2-2.  For 
a given  segment,  the  averaging  process  was  conducted  as  follows. 
During  each  computing  cycle,  the  Lagrangian  nodes  in  contact 
with  the  segment  were  multiplied  by  At,  and  the  sums 

n Fx  At  ’ n FynAtl1  were  bui1t  UP  over  a 9iven  time  interval. 

The  average  is  taken  to  be  the  quotient  of  these  sums  and 
the  time  interval.  The  time-averaged  values  in  each  segment 
after  the  steady-state  was  reached  were  determined,  and  are 
plotted  in  Figures  4.2-1  through  4.2-10. 

The  plots  present  the  time-averaged  force  on  each  sequent, 
the  odd-numbered  plots  giving  the  axial  components  and  the  even- 
numbered  plots  giving  the  radial  component.  The  axial  forces 
are  presented  because  they  represent  the  relative  contribution  of 
each  segment  to  the  projectile's  deceleration.  The  forces  on  the 
conical  nose  tip  are  not  plotted  in  the  Figures  because  they  are 
so  large  that  their  plotting  would  produce  poor  scaling.  Instead, 
they  are  called  out  in  the  figure  captions.  The  area  of  the 

9 

conical  nose  tip  is  7.1659  cm  . 
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To  obtain  the  traction  (force  per  unit  area)  on  a segment, 
the  force  must  be  divided  by  the  segment's  area.  The  areas  are 
tabulated  in  Table  4.2-2,  and  the  tractions  are  presented  in 
Figures  4.2-11  through  4.2-20.  The  tractions  are  the  physical 
entities  which  provide  the  boundary  condition  for  the  stress  field 
in  the  projecti le. 

As  the  soil  flows  around  the  shoulder  of  the  conical  tip, 
it  separates  from  the  projectile  because  of  its  inertia.  This 
creates  a free  surface  on  the  soil  on  which  there  is  no  traction, 
and  the  stress  within  the  soil  causes  the  surface  to  accelerate 
back  to  the  projectile.  This  phenomenon  leads  to  low  average 
traction  on  the  projectile's  surface  just  above  the  conical  tip, 
followed  by  a spike  where  the  soil  reattaches  to  the  vehicle. 

The  region  of  separation  may  be  seen  in  the  configuration  plots 
of  Section  4. 3. 

Because  the  cells  are  given  too  large  a momentum  by  the 
conical  tip  (as  discussed  above),  the  region  of  separation  is 
overestimated.  This  may  be  verified  by  comparing  the  force 
profiles  for  the  coarse  and  fine  zone  cases  shown  in  Figures  4.2-3 
and  4.2-9.  The  coarse  zoning  leads  to  a wider  region  of 
separation,  and  a higher  spike  in  force  (and  therefore  traction) 
when  reattachment  occurs. 

The  forces  and  tractions  which  appear  on  the  conical  nose 
tip  are  inordinately  large.  The  reason  is  that  the  zone  size 
of  all  the  cases  is  too  large  to  resolve  the  true  stress  on  the 
conical  tip.  When  the  tip  first  encounters  the  upper  left  hand 
corner  of  a zone  on  axis,  (See  the  zone  configuration  plots  in 
Figures  4.3-1  through  4.3-5)  that  corner  is  imparted  in  one  or 
two  time  cycles  the  velocity  required  for  it  to  travel  along  the 
tip.  For  an  impact  velocity  of  500  ft/s,  the  increments  in 
x and  y velocity  are  on  the  order  of  600  cm/s.  For  the  purposes 
of  computing  momentum,  each  vertex  in  the  computing  grid  has 
associated  with  it  a quarter  of  the  total  cell  mass  of  each  ad- 
joining cell,  so  that  the  momentum  change  in  the  cell  is 
6000(M/2)  gm  cm/s,  since  there  are  two  adjoining  cells  for  each 
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vertex  on  axis.  The  force  on  the  nose  tip  was  computed  as  the 
time  rate  of  change  of  momentum  it  created  in  the  computing  grid. 
When  the  dimensions  of  a quarter-cell  exceed  the  width  of  the 
conical  nose  tip,  as  they  do  in  all  cases  treated  here,  the  tip 
accelerates  too  much  mass  and  the  stress  and  force  on  the  nose 
tip  is  grossly  overestimated,  because  the  computing  algorithms 
overestimate  the  momentum  increase  for  the  reasons  discussed 
above. 

The  major  point  of  this  discussion  is  that  to  compute  the 
stress  distribution  on  a blunted  vehicle  accurately,  zone 
sizes  of  2 , or  perhaps  even  1 cm  should  be  used.  This  consider- 
ation is  particularly  important  because  the  stress  distribution  is 
required  to  compute  the  structural  response  of  the  EPW. 

The  above  discussion  can  be  made  more  quantitative  by 
examining  the  forces  and  tractions  computed  in  the  various  cases 
of  this  study.  These  are  summarized  in  Table  4.2-1. 


TABLE  4.2-1 


SUMMARY  OF  FORCES  AND 
TRACTIONS  ON  CONICAL  NOSE  TIP 


Case 

Radial 

Force 

Dynes 

Axial 

Force 

Dynes 

C Contribution 
of  nose  tip 
to  Deceleration 

Axial 

Traction 

kbar 

Radial 

Traction 

kbar 

1 

-7.0  x 109 

9.7  x 109 

28.7 

1.35 

-.977 

2 

-6.4  x 109 

6.4  x 109 

40.1 

.888 

-.888 

3 

-6.1  x 109 

6.1  x 109 

58.1 

.852 

-.852 

4 

-3.0  x 1011 

3.0  x IQ11 

83.3 

42.0 

-42.0 

5 

-4.9  x 10° 

4.9  x 109 

33.3 

.678 

-.678 

The  table  shows  that  for  all  the  frictionless  cases 
(2  through  5),  the  radial  and  axial  force  are  the  same,  because  the 
frictionless  45°  tip  imparts  equal  x and  y velocity  increments 
to  the  soil  it  encounters.  Cases  2 and  5 are  different  only  in 
that  Case  5 had  3 cm  zones  and  Case  2 had  4.1275  cm  zones. 

g 

The  more  finely  zoned  case  has  an  axial  force  of  4.9  x 10 

g 

dynes,  compared  with  the  6.4  x 10  dynes  for  coarser  zoning. 

The  reduction  in  force  is  23?'.,  and  is  produced  by  a zone  size 
reduction  of  27%.  The  axial  forces  further  back  on  the  ogive  are 
given  in  Figures  4.2-3  and  4.2-9,  and  it  may  be  seen  that  they 
are  quite  similar  at  distances  greater  than  12  cm  above  the  tip. 

Case  4 shows  that  this  phenomenon  will  have  a very  important 
influence  on  the  computed  deceleration  at  high  velocity,  since 
the  force  on  the  conical  tip  is  responsible  for  83.3%  of  the 
computed  deceleration.  It  seems  likely  that  the  computed 
deceleration  for  this  case  is  seriously  in  error,  and  that  smaller 
zones  should  be  used.  This  question  could  be  resolved  by  using 
successively  finer  zones  until  convergence  to  a common  deceleration 
was  obtained. 


TABLE  4.2-2 


AREA  OF  PROJECTILE 
SEGMENTS  USED  TO  MONITOR  STRESS 


— 

Lower  Station 
of  Segment 
cm 

Segment 

Area 

cm2 

Lower  Station 
of  Segment 
cm 

Segment 

Area 

crrr 

1.27 

12.2493 

24.13 

53.5337 

2.54 

15.4303 

25.40 

54.9265 

3.81 

18.4963 

26.67 

56.2326 

5.08 

21.4496 

27.94 

57.4527 

6.35 

24.2927 

29.21 

58.5878 

7.62 

27.0275 

30.48 

59.6386 

8.89 

29.6562 

31  .75 

60.6057 

10.16 

32.1806 

33.02 

61 .4898 

11.43 

34.6027 

34.29 

62.2913 

12.70 

36.9240 

35.56 

63.0111 

13.97 

39.1464 

36.83 

63.6491 

15.24 

41  .2715 

38.10 

64.2061 

16.51 

43.3005 

39.37 

64.6823 

17.78 

45.2350 

40.64 

65.07801 

19.05 

47.0764 

41  .91 

65.3935 

20.32 

48.8258 

43.18 

65.6290 

21.59 

50.4844 

44.45 

65.7847 

22.86 

52.0534 
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AXIAL  FORCE  (10  DYNES) 


RADIAL  FORCE  (10  DYNES) 


CPU  PRGBLEN  1 


Impact  Speed:  500  ft/s 
Soil  type:  Layer  1 
Friction?:  Yes 
Initial  zone  size:  4.13  cm 
Averaging  interval:  3 - 8 ms 


\ i 


— 10  cm. 


DISTANCE  FROM  NOSE  (10  CM) 


Figure  4.2-2.  Time-averaged  radial  force  profile  for  Case  1. 

Force  on  conical  tip  = 6.978  x lo"  dynes. 
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AXIAL  FORCE  (10  DYNES) 


EPU  PROBLEM  2 


i o. 


DISTANCE  FROM  NOSE  (10  CM) 


Figure  4.2-3.  Time-averaged  axial  force  profileqfor  Case  2. 

Force  on  conical  tip  = 6.365  x lCr  dynes. 
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RADIAL  FORCE  (10  DYNES) 


EPU  PROBLEM  2 


Figure  4.2-4.  Time-averaged  radial  force  profile  for  Case  2. 

Force  on  conical  tip  = 6.365  x 10y  dynes. 
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RADIAL  FORCE  (10  DYNES) 


EPW  PROBLEM  3 


Impact  speed:  500  ft/s 
Soil  type:  Fill 
Friction?:  No 
Initial  zone  size:  4.13  cm 
Averaging  interval:  3 - 0 ms 
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Figure  4.2-6.  Time-averaged  radial  force  profile  for  Case  3. 

Force  on  conical  tip  - 6.103  x IQ9  dynes. 
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AXIAL  FORCE  (10  DYNES) 


EPU  PROBLEM  4 


Impact  Speed:  3500  ft/s 

Soil  type:  Fill  1 

Friction?:  No 

Initial  zone  size:  4.13  cm  j 

Averaqinq  interval:  .5  - 1.16  ms  1 

! I ; 


10  cm. 


EPW  PROBLEM  4 


DISTANCE  FROM  NOSE  UO  ’ CM) 

Figure  4.2-8.  Time-averaged  radial  force  profile  for  Case  4. 

Force  on  conical  tip  = 3.01  x 10^  dynes/ 


AXIAL  FORCE  (10  DYNES) 


EPU  PROBLEM 
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Figure  4.2-9.  Time-averaged  axial  force  profile„for  Case  5. 

Force  on  conical  tip  = 4.856  x lCr  dynes/ 
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RADIAL  FORCE  ( 1 0 DYNES) 


EPW  PROBLEM  5 
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RADIAL  TRACTION  (10  BARS) 


EPU  PROBLEM  1 


Impact  speed:  500  ft/s 
Soil  type:  Layer  1 
Friction?:  Yes 
Initial  zone  size:  4.13  cm 
Averaging  interval:  3 - 0 ms 

: i 


i 


DISTANCE  FROM  NOSE  MO  CM) 

Figure  4.2-11.  Time-averaged  radial  traction  profile  for  Case  1. 
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AXIAL  TRACTION  (10  BARS) 


EPW  PROBLEM  1 


; i Impact  speed:  500  ft/s 

! | Soil  type:  Layer  1 

j : Friction?:  Yes 

Initial  zone  size:  4.13  cm 
i Averaninn  interval:  3 - 3 ms 


B. 


DISTANCE  FROM  NOSE  U0  ' CM) 

Fiqure  4.2-12.  Time-averaoed  axial  traction  profile  for  Case  1. 


EPW  PROBLEM  2 


Impact  speed:  500  ft/s  ! 

Soil  type:  Layer  1 ! 

Friction?:  No 
Initial  zone  size:  4.13  cm 
Averaging  interval:  3 - G ms  ! 
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DISTANCE  FROM  NOSE  U0  ’ CM) 


Figure  4.2-13.  Time-averaged  axial  traction  profile  for  Case  2. 


AXIAL  TRACTION  (10  BARS) 


EPW  PROBLEM  2 


I 


►,  | | Impact  speed:  500  ft/s 

i\  ; ! Soil  tyDe:  Layer  1 

| \ j ; Friction?:  No 

j i ; i Initial  zone  size:  4.13  cm 

i1  ! i i Averaging  interval:  3 - 3 ms 


— :10  cm. 


distance  prom  NOSE  U0  ' CM) 

Figure  4.2-14.  Time-averaged  radial  traction  profile  for  Case  2. 
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RADIAL  TRACTION  (10  ’ BARS) 


EPU  PROBLEM  3 


r " 


ImDact  speed:  500  ft/s 
Soil  type:  Fill 
Friction?:  No 
Initial  zone  size:  4.13  cm 
Averaging  interval:  3 - 0 ms 


DISTANCE  FROM  NOSE  U0  ’ CM) 


Figure  4.2-15.  Time-averaged  axial  traction  profile  for  Case  3. 
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AXIAL  TRACTION  MO  BARS) 


EPU  PROBLEM  3 


Impact  speed:  500  ft/s 
Soil  type:  Fill 
Friction?:  No 
Initial  zone  size:  4.13  cm 
Averaging  interval:  3 - 0 ms 


DISTANCE  FROM  NOSE  MO  ' CM) 


Figure  4.2-16.  Time-averaged  radial  traction  profile  for  Case  3. 
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RADIAL  TRACTION  MO  BARS) 


EPU  PROBLEM  4 


Impact  speed:  3590  ft/s 

Soil  type:  Fill 

Friction?:  No 

Initial  zone  size:  4.13  cm 

Averaging  interval:  .5  - 1.16  ms 
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Figure  4.2-17.  Time-averaged  axial  traction  profile  for  Case  4. 


AXIAL  TRACTION  (10  " BARS) 


EPW  PROBLEM  4 


Impact  speed:  3500  ft/s 

Soil  type:  Fill 

Friction?:  No 

Initial  zone  size:  4.13  cm 

Averaging  interval:  .5  - 1.16  ms 
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Figure  4.2-18.  Time-averaged  radial  traction  profile  for  Case  4. 
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Impact  speed:  500  ft/s 
Soil  tyoe:  Layer  1 
Friction?:  No 
Initial  zone  size:  3.0  cm 
Averaging  interval:  3 - 6 ms 
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Figure  4.2-19.  Time-averaged  axial  traction  profile  for  Case  5. 
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RADIAL  TRACTION  (10  BARS) 


EPW  PROBLEM  5 


Impact  speed:  500  ft/s 
Soil  type:  Layer  1 
Friction?:  No 
Initial  zone  size:  3.0  cm 
Averaging  interval:  3 - 6 ms 


I • \ 

I !; 

i -i 

i 

1 ,i  '• 


— no 


cm. 


. ' '<  • 

DISTANCE  FROM  NOSE  CIO  ’cri) 

Figure  4.2-20.  Time-averaged  radial  traction  profile  for  Case  5. 
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4 . 3 ADDITIONAL  PLOTS  OF  . RESULT S 

In  this  section,  additional  computational  results  are 
presented  in  the  form  of  graphs.  The  graphs  labeled  4.3-1 
pertain  to  Case  1,  etc.  A letter  is  appended  after  the  case 
number  to  identify  the  individual  plot.  The  plots  shown  in  each 
set  are: 

a.  Velocity  versus  time 

b.  Velocity  versus  depth 

c.  Computing  grid  configuration  at  2 ms 

d.  Computing  grid  configuration  at  4 ms 

e.  Computing  grid  configuration  at  6 ms 

f.  Computing  grid  conf iguration  at  8 ms 

Figure  4.3-0  presents  the  initial  computing  grid  confiauration. 

For  Case  4,  the  impact  at  3500  f/s,  the  configurations  are  given 
at  times  when  the  penetration  depth  was  about  1,  2,  3,  and  4 feet. 
Note  that  the  configurations  have  been  expanded  by  a factor  of 
two  laterally  to  provide  better  resolution  of  the  small  zones 
near  the  projecti le. 

For  Case  5,  no  configuration  plots  are  available  for  the 
report.  Instead,  a movie  was  made  of  Case  5 usinq  a continuous 
sequence  of  configuration  plots. 
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.3-la.  Penetrator  velocity  history  for  Case  1. 
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Figure  4.3-ld.  Computing  grid  configuration  at  4 ms  for  Case  1. 
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Figure  4 . 3- 1 e . Computing  grid  configuration  at  6 ms  for  Case  1. 
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Figure  4.3-2c.  Computing  grid  configuration  at  2 ms  for  Case  2 
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Figure  4.3-2d.  Computing  grid  configuration  at  4 ms  for  Case  2 


>4~ 


DEPTH  (10  CM 


DEPTH  MO  CM 


EPU  PROBLEM  2 


CYCLE  -H7 


TIME  •=  3 010-03  SEC 


\ 


\ 


\ : 

\ i 

. : 

, \ • 

: , t'. 

\\ 

s o 

il 

; ’ «. 


. ' "t* 

■ ' >•< 
~ - ■ a 


/ 


*• 

>/ 


t* 


\ r 


. l . 


RADIUS  (10  CM.  ) 


Figure  4.3-2f.  Computing  grid  configuration  at  8 ms  for  Case  2. 
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Figure  4.3-3a.  Penetrator  velocity  history  for  Case  3. 
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Figure  4.3-3d.  Computing  grid  configuration  at  4 ms  for  Case  3. 
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Figure  4.3-3e.  Computing  grid  configuration  at  6 ms  for  Case  3. 
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Figure  4.3-3f.  Computing  grid  configuration  at  8 ms  for  Case  3. 
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Figure  4.3-4a.  Penetrator  velocity  history  for  Case  4 
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Figure  4.3-4c.  Computing  grid  configuration  at  2 ms  for  Case  4 
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Figure  4.3-4f.  Computing  grid  configuration  at  8 ms  for  Case  4. 
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Figure  4.3-5a.  Penetrator  velocity  history  for  Case  5. 
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Figure  4.3-5b.  Velocity  penetrator  depth  cross  plot  for  Case  5. 
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5.  COMPARISON  WITH  EXPERIMENTAL  DATA  AND 
OTHER  COMPUTATIONS 


For  the  purposes  of  computing  the  penetrator's  structural 
response  and  deceleration,  the  most  useful  information  which  is 
obtained  from  these  computations  of  rigid  body  penetration  is  the 
surface  tractions  which  appear  on  the  penetrator.  These  comprise 
the  proper  boundary  conditions  for  other  computations  which  seek  to 
compute  the  stresses  and  deformation  which  arise  within  the  vehicle, 
the  wave  propagation  within  it,  and  its  structural  response.  The 
axial  and  radial  components  of  the  traction  are  presented  in  Figures 
4.2-11  through  4.2-20. 

There  is  no  experimental  data  against  which  to  compare  these 
computed  tractions.  The  measurements  of  soil  stress  made  in  the 
firings  were  taken  so  far  from  the  vehicle  that  they  are  very 
insensitive  to  the  stresses  on  the  vehicle  itself.  The  stress 
field  in  the  soil  decreases  extremely  rapidly  with  radial  distance 
from  the  vehicle.  The  soil  stress  measurements  were  taken  at  a 
stress  level  of  less  than  3 bars  in  layer  1,  and  it  appears  hopeless 
to  attempt  to  use  them  as  a base  from  which  to  extrapolate  back  + 
the  vehicle  to  find  the  forces  on  it. 

The  only  other  experimental  data  is  the  velocity-time  history 
of  the  vehicle,  which  provides  a measure  of  its  deceleration.  Since 
the  deceleration  is  given  by  the  integrated  axial  force  on  the  vehicle 
divided  by  its  mass,  comparing  the  computed  and  measured  deceleration 
has  some  relevance  to  the  computed  stresses.  This  comparison  is  conducted 
in  Section  5.1.  In  Section  5.2  the  computed  tractions  are  compared 
with  those  resulting  from  previous  computations. 
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5.1  COMPARISON  OF  THE  MEASURED  AND  CALCULATED  DECELERATION 


The  deceleration  was  measured  by  accelerometers  which  were 
carried  on  board  the  penetrator.  The  data  were  read  out  in  two  ways. 
First,  the  data  was  transmitted  out  in  real  time  during  the  penetration 
process.  Second,  the  data  was  recorded  in  an  onboard  memory  unit 
during  penetration,  and  was  telemetered  out  after  the  penetrator  came 
to  rest.  The  two  sets  of  data  contain  oscillations  which  are  presently 
unexplained. 

n 3i 

As  our  basis  for  comparison,  we  use  the  smoothed  set  of  dataL  J 
which  is  presented  in  Figure  5.1-1.  In  that  figure,  the  data  is 
approximated  by  straight-line  segments  which  give  the  best  fit  to  the 
onboard  memory  data  which  has  been  filtered  to  200  Hz. 

The  decelerations  computed  in  this  study  are  61.2  g's  for  the 
fill  material  (Case  3)  and  81.9  g’s  for  the  material  of  layer  1 
(Case  5).  WES  has  identified  that  the  top  3 to  5 feet  of  the  firing 
site  are  composed  of  fill  material,  with  the  most  representative  depth 
being  4 feet.  The  next  few  feet  beneath  it  are  of  layer  1 material. 
These  decelerations  are  plotted  in  Figure  5.1-1  as  the  heavy  dashed 
lines,  and  are  to  be  compared  with  the  measured  deceleration  for 
firing  06,  which  is  shown  as  the  heavy  solid  line  in  the  figure.  The 
computed  decelerations  compare  favorably  with  the  measured  values,  but 
are  somewhat  higher  for  both  media. 

As  discussed  in  detail  in  Section  4.1,  zone  sizes  of  3 or  4 
centimeters  are  too  large  to  resolve  the  flow  about  the  blunted  tip, 
and  their  use  leads  to  values  of  deceleration  which  are  too  high. 

To  determine  the  true  decelerations,  the  calculations  should  be  rerun 
using  smaller  zones.  If  this  were  done,  many  useful  results  would  be 
obtained.  If  the  correctly-computed  decelerations  proved  to  agree 
precisely  with  the  measured  values,  it  would  indicate  that  the  material 
property  measurements  conducted  in  WES  provided  sufficient  data  to 
compute  penetrator  decelerations  in  these  soil  types.  In  particular, 
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Impact 


Comparison  between  the  computed  and  measured  decelerations  for  firing  06  at  the 
ORES.  (Note  that  the  measured  deceleration  exhibits  a general  increase  with 
depth  and  peaks  near  12  feet.  The  two  layers  computed  reflect  this  increase, 
but  the  material  near  12  feet  must  be  used  in  a computation  to  obtain  a comparison 
with  the  measured  peak. ) 
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it  would  not  be  necessary  to  include  the  effects  of  friction  in  the 
calculation.  If,  on  the  other  hand,  the  correctly-computed  decelerations 
turned  out  to  be  less  than  the  measured  values,  one  could  determine 
the  amount  of  friction  which  has  to  be  included.  It  is  recognized,  of 
course,  that  such  conclusions  assume  that  the  decelerations  were 
accurately  measured.  If  this  turns  out  not  to  be  true,  the  conclusions 
would  be  deferred  until  correct  measurements  were  taken. 

Before  the  present  study  was  conducted,  the  opinion  had  been 
expressed  by  others  that  it  was  not  economically  feasible  to  utilize 
zones  much  smaller  than  3 cm.  It  was  feared  that  using  small  zones 
would  require  prohibitively  large  amounts  of  computer  time.  However, 
our  research  has  resulted  in  greatly  increased  computing  efficiency. 

For  example.  Case  5,  in  which  three  feet  of  penetration  were  computed 
using  3 cm  zones,  required  less  than  $80  of  computer  time.  Using 
smaller  zones  to  compute  the  correct  deceleration  is  therefore  eminently 
feasible. 

The  measured  decelerations  in  Figure  5.1-1  show  that  a peak 
of  deceleration  occurs  in  the  interval  between  8 and  13  feet.  Peaks 
at  about  this  same  depth  are  present  to  some  degree  in  all  four  of  the 
tests  conducted  in  this  area,  namely.  Tests  01,  02,  06,  and  07,  whose 
parameters  are  called  out  in  the  legend  of  Figure  5.1-1.  Data  taken 
from  nearby  borings,  and  the  pit  which  was  excavated  postshot, ^>14] 
show  that  a coarse-grained  sandy  layer  occurs  near  a depth  of  12  feet. 

It  would  be  useful  to  conduct  a computation  using  the  properties  of  that 
layer,  to  see  whether  the  deceleration  peak  would  be  reproduced  by  the 
computation.  If  it  were,  it  would  add  further  weight  to  our  hypothesis 
that  our  computation,  using  laboratory  data  as  an  input,  can  predict 
the  deceleration.  The  calculations  in  the  first  two  layers  show  an 
increasing  deceleration  with  depth.  A successful  computation  of  the 
deceleration  near  12  feet  would  complete  a tentative  explanation  of 
the  increase  of  measured  deceleration  with  depth,  peaking  at  about 
12  feet. 
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5.2  COMPARISONS  WITH  OTHER  CALCULATIONS 


n 4i 

In  reviewing  the  previous  computations,  J it  was  found  that 
no  two  computations  had  been  conducted  with  the  same  inputs.  The 
most  similar  computations  had  been  conducted  by  SLA  and  CRT.  However, 
CRT  had  used  a blunted  tip,  a rigid  body  approximation,  and  an 
artificially-constructed  model  of  friction.  SLA  had  conducted  two 
usable  computations.  The  first,  which  we  call  SLA-R,  used  a rigid 
body,  but  had  assumed  a sharp-tipped  complete  ogive  configuration  for 
the  nose.  Their  second  calculation,  called  SLA-D,  used  the  blunt-tipped 
configuration,  but  assumed  that  the  body  was  deformable.  Other 
differences  were: 

1.  Both  SLA  cases  used  a different  fit  of  the  soil 
properties  than  CRT,  and  the  fits  reflected  substantially 
different  soil  properties. 

2.  Neither  SLA  computation  included  friction. 

These  inputs  are  summarized  in  Table  5-1. 


TABLE  5-1 

INPUTS  TO  COMPUTATIONS 

(All  these  computations  used  Lagrangian  codes) 


CRT 

SLA-R 

SLA-D 

Case  1 

Case  5 

Rigid  Body 

yes 

yes 

no 

yes 

yes 

Blunt  Tip 

yes 

no 

yes 

yes 

yes 

Friction 

yes 

no 

nc 

yes 

no 

Zone  Size 

4 cm 

3 cm 

3 cm 

4 cm 

3 cm 

Material  Model 

CRT  1 

SLA  1 

SLA  1 

CRT  1 

CRT  1 

CRT  1 denotes  the  CRT  fit  to  layer  1 material. 
SLA  1 denotes  the  SLA  fit  to  layer  1 material. 
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For  these  reasons  the  computational  review1-  J concluded  that 
it  was  not  relevant  to  conduct  detailed  comparisons  between  the  CRT 
and  SLA  computations. 

Case  1 of  the  present  study  was  designed  to  simulate  the  CRT 
calculation  as  closely  as  possible.  It  is  therefore  relevant  to 
conduct  comparisons  between  the  surface  tractions  calculated  in  Case  1 
and  by  CRT.  This  is  done  in  Section  5.2.1. 

Case  5 of  this  study  uses  the  same  vehicle  configuration  as 
SLA-D,  and  the  same  zone  size.  Like  the  SLA  case,  no  friction  is 
included.  The  major  differences  between  the  two  computations  are  that 
SLA-D  used  a deformable  body,  and  a different  soil  model.  Therefore, 
it  is  less  relevant  to  compare  Case  5 with  SLA-D,  but  some  qualitative 
comparisons  are  conducted  in  Section  5.2.2. 


5.2.1  Comparison  with  the  CRT  Calculation 


The  inputs  of  Case  1 simulated  those  of  the  CRT  calculation  as 
closely  as  possible.  It  is  therefore  relevant  to  conduct  comparisons 
between  the  stresses  calculated  in  the  two  studies. 

Both  calculations  used  the  same  vehicle  geometry  (blunted-tip) , 
same  soil  material  (in  fact  the  CRT  fit  was  used  in  Case  1),  the 
same  zone  size,  and  the  same  physical  model  of  friction.  Specifically, 
the  shearing  stress  on  the  vehicle  is  set  to 
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Figure  5.2-1  compares  the  normal  stress*  on  the  vehicle  given 
by  the  two  computations,  and  Figure  5.2-2  compares  the  computations  of 
tangential  stress.*  The  results  for  the  present  study  were  obtained 
by  summing  the  appropriate  components  of  the  axial  and  radial  tractions 
shown  in  Figure  4.2-11  and  12. 
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♦While  the  terms  "normal  stress"  and  tangential  stress"  are  coming  into 
corrmon  usage,  they  actually  refer  to  the  normal  and  tangential  components 
of  surface  traction. 
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To  1640  bars  at 
y = .635  cm 


normal  component  of  surface  traction  as  computed  in  this 
and  by  CRT.  [10] 


to  269  bars  at  y = .635  cm 


Axial  Distance  from  Tip  (cm) 

son  of  tangential  component  of  surface  traction  as  computed  in  this 
Case  1)  and  by  CRT.  [10] 


It  may  be  seen  that  the  results  are  different  in  several  respects. 
First,  the  normal  stress  computed  in  this  study  is  in  general  higher 
than  that  computed  by  CRT.  The  most  important  reason  for  this  is  that 
the  present  study  utilizes  an  algorithm  which  provides  a stress  on  the 
vehicle  itself,  whereas  the  CRT  computation  used  the  stresses  in  the 
centers  of  the  computing  cells  adjacent  to  the  vehicle.  Because  the 
stress  in  the  soil  varies  extremely  rapidly  with  distance  from  the  vehicle, 
this  tends  to  cause  CRT  to  underestimate  the  stress.  Secondly,  the  computing 
algorithms  utilized  by  CRT  did  not  permit  the  soil  to  separate  from  the 
vehicle  as  it  flowed  around  the  blunted  tip.  Thus,  the  region  of  small 
stress  which  occurs  just  above  the  conical  tip  in  the  present  study 
is  not  reflected  in  the  CRT  result. 

The  agreement  between  the  two  computations  of  stress  becomes 
better  with  increasing  distance  from  the  penetrator  tip.  The  agreement 
becomes  very  good  at  about  30  cm  from  the  tip,  and  both  computations 
show  that  little,  if  any,  stress  is  present  at  distances  greater  than 
40  cm  from  the  tip.  Inasmuch  as  the  stress  gradients  in  the  soil  are  very 
large  near  the  tip,  and  become  smaller  as  one  moves  back  along  the 
penetrator,  this  variation  in  agreement  is  precisely  what  one  would  expect. 

In  the  CRT  computation,  computing  nodes  which  encounter  the  blunted 
tip  are  imparted  the  velocity  necessary  to  keep  them  on  the  tip.  This 
procedure  increases  the  momentum  in  the  soil,  but  does  not  decrease 
the  momentum  in  the  penetrator.  The  decrease  in  penetrator  momentum 
comes  only  when  the  motions  produced  in  the  soil  generate  strains, 
which  results  in  stresses  in  the  computing  grid.  The  net  effect  is  to 
produce  more  axial  momentum  in  the  soil  than  the  penetrator  has  lost. 

The  final  result  is  that  the  decelerate  of  and  the  tractions  on  the 
penetrator  are  underestimated. 

One  would  therefore  expect  that  reducing  the  grid  size  in  the 
CRT  computations  would  lead  to  larger  computed  decelerations  and  larger 
tractions  on  the  penetrator,  because  the  errors  which  cause  the  under- 
estimate are  reduced  thereby.  For  the  reasons  discussed  in  Section  4.2, 
reducing  the  grid  size  in  the  present  study  would  lead  to  smaller 


decelerations,  smaller  tractions  on  the  conical  tip,  and  a different 
traction  distribution  on  the  ogive.  It  is  therefore  anticipated  that 
the  decelerations  computed  by  CRT  and  those  of  the  present  study  provide 
upper  and  lower  bounds  on  the  correct  deceleration. 

5.2.2  Comparison  with  the  SLA  Computation 

In  this  section,  we  discuss  SLA  deformable  body  calculation, 
which  is  Case  SLA-D  in  Table  5.1.  The  inputs  for  this  computation  were 
discussed  in  Section  5.2. 

In  the  SLA  deformable  body  calculation,  the  interaction  between 
the  soil  and  the  penetrator  is  accomplished  using  a slide  line  formulation 
similar  to  that  of  Wilkins.  The  slide  line  formulation  permits 
differential  motion  between  the  soil  and  penetrator,  but  computes  the 
stresses  which  arise  in  both  media.  The  formulatic.  used  by  SLA  does 
not  permit  the  soil  to  separate  from  the  penetrator  Therefore,  as  in 
the  CRT  calculation,  soil  nodes  which  encounter  the  proiectile  are  imparted 
the  velocity  necessary  to  keep  them  on  the  blunted  tip.  The  essential 
difference  is  that  the  stress  in  the  penetrator  required  to  produce  the 
soil  acceleration  is  taken  into  account.  The  stresses  are  the  agency 
which  impart  momentum  to  the  soil,  and  by  Newton's  third  law  impart  an 
equal  and  opposite  momentum  to  the  projectile.  Therefore,  in  the 
SLA  computation,  the  momentum  acquired  by  the  soil  comes  at  the  expense 
of  a momentum  loss  of  the  projectile.  The  computing  algorithms  of  the 
calculations  presented  here  accomplish  the  same  end,  even  though  a rigid 
body  is  used. 

The  SLA  computations  used  3 cm  computing  cells.  For  the  reasons 
discussed  in  Section  4.2,  use  of  these  large  computing  cells  causes  the 
momentum  imparted  to  the  soil  to  be  overestimated,  which  tends  to 
overestimate  the  vehicle  deceleration.  In  addition  the  large  cells  should 
generate  stresses  on  the  blunted  tip  of  the  same  order  as  those  which 
were  computed  in  the  present  study.  As  shown  in  Table  4.2-1,  the  time- 
average  value  of  the  surface  axial  surface  traction  for  Case  5 is 
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.678  kbars,  and  substantial  oscillations  about  the  average  value 
occurred.  The  oscillations  in  traction  are  generated  as  each  cell 
encounters  the  blunted  tip,  and  we  would  expect  the  same  phenomena  to 
be  present  in  the  SLA  computation.  This  would  cause  large  oscillation 
of  the  stress  and  strain  fields  within  the  projectile.  An  examination 
of  the  computed  SLA  data  confirms  that  this  is  true.  The  oscillations 
in  the  velocity  field  within  the  penetrator  were  sufficiently  large  to  mask 
the  computed  deceleration.  In  addition,  mean  stresses  as  high  as  1.02 
kbars  and  octahedral  shear  stresses  as  high  as  2.9  kbars  were  computed 
within  the  penetrator. 

Because  only  the  pressure  (mean  stress)  and  octahedral  shear 
stress  are  available  in  the  SLA  results,  it  is  not  possible  to  determine 
the  surface  tractions  and  conduct  a comparison. 

Our  interpretation  of  these  results  is  that  much  finer  zoning 
in  the  soil  would  be  required  in  the  SLA  computations  to  compute  the 
correct  stresses  within  the  vehicle  and  its  deceleration. 
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